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ABSTRACT 

This study is an investigation of the stellar density profile of the Galactic disc in the Anti- 
centre direction. We select over 40,000 early A stars from IPHAS photometry in the Galac- 
tic longitude range 160° < / < 200° close to the equatorial plane (-1° < b < +1°). We 
then compare their observed reddening-corrected apparent magnitude distribution with sim- 
ulated photometry obtained from parameterised models in order to set constraints on the 
Anticentre stellar density profile. By selecting A stars, we are appraising the properties 
of a population only ^100 Myrs old. We find the stellar density profile of young stars 
is well fit to an exponential with length scale of (3020 ± l2Q stat i S ti C ai ± systematic) pc, 
which is comparable to that obtained in earlier studies, out to a Galactocentric radius of 
Rt = (13.0 ±0.5 statistical ±0.6 'systematic) kpc. At larger radii the rate of decline appears to 
increase with the scale length dropping to (1200 ± 30Q s tatisticai ±" '0 systematic) P c - This result 
amounts to a refinement of the conclusions reached in previous studies that the stellar den- 
sity profile is abruptly truncated. The IPHAS A star data are not compatible with models that 
propose a sudden change in metallicity at Rq = 10 kpc. 

Key words: surveys - Galaxy: abundances - Galaxy: disc - Galaxy: structure - Galaxy: 
stellar content 



1 INTRODUCTION 

Our position within the Milky Way Galaxy makes studying its 
structure and evolution both uniquely challenging and rewarding. 
Its constituent components can be observed in unique detail. How- 
ever, disentangling its structure is complicated by difficulties in ob- 
taining accurate distances to objects and frequently heavy extinc- 
tion. Studying the Galactic thin disc is particularly difficult; within 
the thin disc extinction is at its strongest and its constituent stellar 
population is a smooth mix of stars formed over a wide range of 
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ages, in vastly differing conditions and locations, and which mi- 
grate across the Galaxy. 

Such difficulties help to explain the vast range of re- 
sults which have been obtained when determining parameters 
such as the scale length of the stellar density in the thin 
disc. Early estimates tended to be relatively long, for example: 
Ide Vaucouleurs & Pence 1 1978 ) obtained 3800 pc, van der Kruit 
( Il986h 5500 pc, lHabind dl988h 4500 pc. iLewis & Freemar 
dl989h 4400 pc and kent. Dame & Faziol dl99ll) 3.0 kpc. 
iRobin. Creze & Mohanl d 19921) compared their Galactic model to 
CCD photometry of a 29 arcmin 2 area in SA23 and found a 
considerably shorter scale length of 2500 pc. Ruphv et al. ( 1996) 
adopted a similar approach, but this time using photometry from 
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the DENIS survey , 
Robi n et al. I d2003h obta i ned a scale length of 2530 pc from DE 
NIS data. Freudenreich d 19961) derived a length of 2600 pc from 
DIB R E observat i ons.lBuser. Rons & Karaalil d 19981) found 4000 pc 
an d lSiegel et alJd2002h~ found 2500-3125 pc, once a correcti on for 
binarity has been applied. Most recently [juric et all d2008h mea- 
sured the thin disc scale length, using high Galactic latitude pho- 
tometry from SDSS, finding a scale length of 2600 pc. Many of the 
previous studies have not used observations at low latitude, reduc- 
ing the sensitivity of their measurements and preventing them from 
studying the disc at larger Galactocentric radii. 

It is important to note that the scale length of the thin 
disc in spiral galaxies is expected to change: it may increase 
with time due to so called inside-out galaxy formation. This 
process involves the galaxy accreting gas in its outer disc at 
later times, from which it forms new stars, whilst simultane- 
ously exhausting gas supplies i n its core. Inside-out forma - 
tion is discussed inlLarsoa dl976h. iMatteucci & Francois! dl989h . 
Chiap pini. Matteucci & Grattonl J 1997b and many others. Within 
the local group, it has been suggested that inside-out formation 
would accou nt for observed metallicity and mean stellar age gradi- 
ents in M33 (Ma grini et ail2 007; Verlev et al. 2009; Williams et al 



2009 ) . In our own Galaxy, several authors (inclu ding: ICescutti et al 



20071; IColavitti et al.l2009l ; iMagrini et al]|2009h have demonstrated 



that chemodynamical models of the formation and evolution of the 
Galaxy, which incorporate inside-out formation, produce results 
which are consistent with some observations of radial abundance 
gradients. 

The effect of studying a Galaxy with a thin disc scale length 
which varies with time is that any measurement of the scale length 
from stars will be dependent on the types of stars used. If stars could 
be easily grouped by age, one would expect the younger groups 
to show a longer scale length than the older ones if the Galaxy 
was undergoing inside-out formation. Unfortunately, in practice it 
is not possible to determine accurate ages for most stars and so 
any samples selected on the basis of spectral type for example will 
contain a superposition of ages and thus the measured scale length 
will be dependent on the range of ages covered by the sample. 

Ivan der Kruitl d 19791) first observed that the exponential discs 
of some spiral galaxies are radially truncated, Indeed it has subse- 
quently been observed that truncation is present in the light pro- 
files of the majority of spiral galaxies jPohlen. Liitticke & DettmaJ 
l200ll ; iKregel, van der Kruit & de Griislp2002T) . With local galaxies 
it is possible to resolve individual stars an d so study the radial 
distribution of a giv en range of stars. M33 (Ferguso n et al.ll2007l ; 
Williams et al . 2009) shows truncation in the radial distribution of 
predominantly RGB stars. Other local galax i es, such as NGC 300 
iBland-Hawthorn. Vlaiic. Freeman & Drainej d2005h . appear not to 
be trunca ted. In this context , the truncation of the G alactic disc, as 
shown bv lRobin et aT] dl992h and lRuphv et"al] J 19961) is not surpris- 
ing. The key difference between their results and studies of other 
galaxies, is that they find that their observations are well described 
by a model Galactic disc which is sharply truncated, with stellar 
density dropping to zero beyond some Galactocentric radius, rather 
than a model where the truncation is less sharp. It is worth noting, 
that although a model with such a sharp truncation is a good fit to 
their observations and a better fit than no truncation, it is possible 
that a mod el with a les s sharp truncation would also fit their ob- 
servations. iRuphv et alj dl99efc find that the truncation occurs at a 
radius of 15 ± 2 kpc, whilst iRobin et al] d 19921) find that a model 
with truncation at a radius of 14 kpc best describes their observa- 
tions. 



The caus e of radial truncation in gal actic discs is poorly un- 
derstood (e.g. lElmegreen & Hunted [2006 ) . There have been sev- 
eral mechanisms proposed to explain it. IKennicuttl d 19891) sug- 
gested that star formation in a gaseous disc may be truncated 
when the column density of the gas drops bellow the threshold 
for gravitational instabilities to form, based on t he sharp trun- 
cation of Ha emission observed in many discs. iBoissier et al.l 
d2Q03[) were able to show that truncation of the form proposed 
by IKennicuttl d 1989b was consistent with their observations of the 
Galaxy. iNaab & OstrikeJ (2006) suggested that truncation due to 
this mechanis m would occur at a Galactocen tric radius of ~ Ylkpc 
in the Galaxy. lElmegreen & Parravanol d 19941) andlSchave (2004) 
suggest that the transition of the ISM from being dominantly in 
the cold phase to being dominantly in the warm phase is the more 
correct physical factor s uppressing star form ation. A very different 
view has been taken by Ivan der Kruitl dl987l) . who suggested that 
the gaseous disc of a galaxy may be truncated during the galaxy's 
initial formation. Furt her mechanisms for truncating the disc in- 
volve magnetic fields (Battaner, Florido & Jimenez-Vicentejl2002l) 
and the outer Linblad resonance dErwin. Pohlen & Beckmanl2008l) 

iRoskar et al.l (2008) demonstrate how the sharpness of disc 
truncation may be reduced by the redistribution of stars as a re- 
sult of secular processes in the disc. Azzolli ni. Truiillo & Beckmard 

(2008) show, using fi-band photometry, that the truncation radii 
of galaxies incre ase with age, whilst the surface brig htness at the 
break decreases. IPflamm-Altenburg & Kroupa d2008l) suggest it is 
possible for a galaxy to exhibit sharp disc truncation in Ha emis- 
sion without this suppression of star formation. 

It is important to be aware of what is viewed as truncated in 
each study. Observations of the Galaxy and other local galaxies are 
able to study the truncation in the numbers of certain types of stars. 
In more distant galaxies the ability to resolve individual stars is lost. 
As a result, it is the light profile of more distant galaxies that is 
freque ntly examined for truncations. BothlMartmez-Serrano et all 

(2009) and lSanchez-Blazquez et al.l d2009l) demonstrate models of 
galaxy formation which can produce galaxies with truncated light 
profiles and untruncated stellar mass profiles. The difference in the 
profiles is the result of the mean age of stars varying with radius, 
beyond the observed truncation radius in the light profiles the mean 
age of stars increases resulting in a smaller number of young and 
bright earl y-type stars. These models are in line with the obser- 
vations of lBakos. Truiillo & Pohlenl d2008l) who observe consistent 
(g' — r') profiles. Therefore, it is necessary to study truncation in the 
discs of galaxies with resolved stellar populations so that the effects 
of changing stellar age and density can properly be understood. Of 
the galaxies in which we can resolve their stellar populations it is 
our own which can be studied in the greatest detail. 

The radial dependence of galactic metallicity is al so unclear. 
Meta ll icity may drop with a constant gradient (e.g. iFriel et all 
120021 ; IBragaglia & Tosil 120061 ; iLuck et alj |2006|). or the gradi- 
ent may flatten at large Galactocentric ra dii dCarraro et all 2007; 
IBragaglia et alj2008l;lLemasle et al.l2008l) . The latter option would 
agr ee with the chemical e volution models of Chiappin i et alj d 19971) 
and lCescutti et alj d2007l) . These studies assume the inside-out for- 
mation of the galactic disc proceeded via two episodes of sig- 
nificant material infal l, the first being assoc iated with the for- 
mation of the Galaxy. Maciel & Costa d2009l) estimated the time 
evolution of the Galactic metallicity gradient, based on central 
stars of planetary nebulae with a range of ages. They conclude 
that the gradient appears to have been steeper in the past, a 
finding which would agree with inside-out formation. An al- 
ternative scheme is that the metallicity obeys a step function, 
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dropping suddenly at Rg ~ 10 kpc I Corder & Twarod l200ll ; 
lYong. Carney & Teixera de Almeidall2005l) . this could arise as the 
result of a significant accretion event. The tracers used to study 
metallicity will affect the result s obtained: the use of re latively old 
objects (e.g. old open clusters in Bragaglia & Tosi 2006) compared 
to younger ones (e.g. OB stars in lDaflon & Cunha 2004) will probe 
the metallicity in the ISM at an earlier age. Thus, as the metallicity 
of the ISM will have evolved over time, the results obtained will 
be different. Even where the choice of scheme agrees and the same 
tracers are us ed there is still signi ficant disagreement on the size of 
any gradient jLemasle et a"il2 008. and references therein). 

Studies of other spiral galaxies have revealed that they too ex- 
hibi t metallicity gradien ts. M3 1 appears to show a metallicity gradi- 
ent dWorthev et al.l2005h . though analysing it is complicated by the 
inclination of M3 1 and difficulties in separating contributions from 
its disc and bulge components. Due to its more favourable inclina- 
tion M33 is frequently considered to be better suited for studying 
abundance g radients and has long been known to posses them (e.g. 
ISearld f 1971). More recently these have been studied in greater de- 
tail, using a range of tracers, for example: planetary nebulae (e.g. 
Magrini et al. 2004); red giant branch stars ^Barker et al.l2007l) : HII 
regions dMagrini et al.ll2007h . As with our Galaxy, the results ob- 
tained f or M33 can diverge , even if the same tracer has been em- 
ployed dMagrini et alj|2007i and references therein). iMagrini et al.l 
(2009) determine that the metallicity gradients in M33 do not ap- 
pear to have changed significantly with time. Outside of the Lo- 
cal Group, metallicity gradients have been recor ded in other spi- 
ral galaxies, such as NGC 300 dVlaiic et al.ll2009h - The possible 
flattening of the metallicity gradient at large radii does not ap- 
pear to be a unique feature of our Galaxy: s imilar behaviour has 
been observed in M31, M83 and NGC 300 dWorthev et aTll2005l ; 
iBresolin et alj|2009l : lvla"uc et alj2009T) . 

We investigate the outer disc by comparing the photometry of 
A stars derived from IPHAS dDrew et al.ll2005l) . with simulations 
exploiting aspects of the Besancon formulation of Galactic models 
dRobin et all2003l) . In section[2]we introduce the CCD photometry 
we have used and detail how it has been prepared. Section [3] ex- 
plains why and how observed and simulated photometry are com- 
pared. The results obtained are presented in section [4] along with 
discussion and analysis of other relevant parameters. 



2 OBSERVATIONS 

The INTAVFC Photometric Ha Survey of the northern Galactic 
Plane (IPHAS; brew et al.ll2005h is the first comprehensive digital 
survey of the disc of the Galaxy (\b\ ^ 5°), north of the celestial 
equator. Imaging is performed in the r', (' and Ha bands down to 
r' ~ 20 (10c). At the time of writing all fields in the survey area 
have been observed at least once, with a minority of these in need 
of replacement by data obtained in better observing conditions. 

In this study we use only observations near the Galactic An- 
ticentre, of fields with their centres in the region 160 < / < 200, 
— 1 < b < 1. Along sightlines towards the Galactic Anticentre 
Galactocentric radius increases at its greatest rate with heliocentric 
distance. As such, stars in the outer Galaxy are brightest along these 
sightlines. T he Galactic Antice ntre is a known node in the Galac- 
tic warp (e.g. lRevle et al.l2009l) . so by choosing sightlines near this 
direction the influence of the warp on observations is reduced. 

It is possible to use the spatial distribution of a particular 
type or types of object to trace Galactic structure. Studies of the 
gaseous component of the disc have utilised the 21cm HI emis- 



sio n line (e.g Kalberla et al. 2005). T o study the stellar density pro- 
file lRobin et alJdl992l) , lRuphv et al.l dl996h and lRevle etai] d2009l) 
have used photometry of all stars. This approach maximises the 
number of stellar tracers available, reducing statistical uncertain- 
ties. T his was necessary give n the relatively sma ll samples avail- 
able to lRobin et ai1dl992h and lRuphv et al.ldl996h . 

The cost of using all stars is that it conceals any distribution 
differences which occur between different types of star. Such con- 
trasts could arise from changes in the structure and composition of 
the outer Galactic thin disc as a function of time. Even if a rela- 
tively restricted subset of stars is studied, a wide range of ages may 
be admitted and so results would be blurred. This is clearly true of 
subsets of late type-stars, which may cover an age range of up to 
10 Gyrs. 

In this study the tracers used are early- A type dwarfs. By se- 
lecting these the age range is minimised, as they only remain on 
the main sequence for ~ 200 Myrs. Therefore, we are measuring 
the distribution of stars that are relatively homogeneous with re- 
spect to age and so can avoid the problems discussed above for 
samples which span a broad range of ages. Early A-stars stars are 
also bright (M r » ~ 1.5) and so the distance range over which they 
can be studied is large. 

The IPHAS (r' — Ha) colou r is strongly corre lated with the 
equivalent width of the Ha line dDrew et ai] |2005l) . Ha absorp- 
tion is strongest in stars with effective temperatures of ~ 9000 K, 
which for solar metallicities corresponds to a spectral type of A2- 
A3. The stronger gravitational broadening of the main sequence 
stars, with respect to evolved stars, further strengthens the Ha 
equivalent width and so separates these classes of stars on an 
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Ha) colour-colour diagram. This separation is not 



achieved with a purely broadband filter set. The effect of extinc- 
tion on the IPHAS colour-colour diagram is such that there is good 
discri mination between reddening and in trinsic (i.e. unreddened) 
colour dDrew et al.ll2"005l;ISale et alj|2009l) . Therefore, as described 
in detail by Dre w et al.l d2008f) . early A-stars, on or near the main 
sequence, are easily identifiable in IPHAS colour-colour plots and 
it is straightforward to measure their ex tinction. 

We adopt the colour cuts of iDrew et al] {2008), with 
8(r' — Ha) = 0.03. Solar metallicity stars selected in this man- 
ner will have spectral types in the approximate range A0-A5. 
Selected stars with super or sub-solar metallicities will formally 
match slightly different spectral types, but their range of intrinsic 
(/ — (') colour and effective temperature will be unaffected. As 
only main sequence stars are included in the sample, the age range 
is relatively narrow, from roughly 10 Myr up to ~ 200 Myr. 

As observations approach the faint magnitude limit of the sur- 
vey, the sample chosen will become progressively more mixed with 
other types of stars as a result of increasing photometric errors. Our 
simulations indicate, that, as expected, most of these contaminat- 
ing stars are later A-type dwarfs. A magnitude limit of r' ^ 19 was 
imposed on the sample, to keep the contamination level integrated 
across the entire sample below ~ 10%. Given an asymptotic Galac- 
tic extinction of A r i = 2.5, typical for low latitude sightlines near 
the Anticentre, and an absolute magni tude of M/ = 1.5, typical for 
stars in the sample dHouk et al.|[l997l) . this magnitude limit corre- 
sponds to a heliocentric distance range of 10 kpc. We found that 
altering the magnitude limit to r < 19.5 did not significantly affect 
the results obtained, despite increasing the overall contamination to 
~ 20%. 

Initially, photometry obtained in observing conditions not 
meeting the survey threshold was removed from the prospective 
sample. Observations were rejected if they exhibited seeing > 1.7", 
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Figure 1. Black points show data from IPHAS field 4199, where 13 ^ r 1 < 
20. Main sequences where extinctions equivalent to Ay = 0, 4 , 8 for an 
AOV star have been applied are shown with solid blue lines and the dashed 
green line shows an A3V reddening line. Stars which would be selected by 
the colour-cut described above are shown with larger red points. 

ellipticity > 0.2 or a lOo limiting magnitude brighter than r' — 20, 
which is indicative of low sky transparency. The seeing and ellip- 
ticity values used are those measure d during the pipeline reductio n 
of IPHAS data (for more details see Gonza lez-Solares et alfeooih . 
This cleaned sample contained observations of 271 fields. 

As the IPHAS survey has yet to be placed on a uniform pho- 
tometric scale there is a minority of IPHAS observations which 
exhibit uncorrected photometric offsets. In such observations the 
loci of all stars, in cluding early A-sta rs, on the measured colour- 
colour plane shifts. iDrew et alj (2005) and Fig. Q] demonstrate that 
early A-stars always lie at the bottom of the main stellar locus. 
Therefore, their locus is easily observed and any offset from the 
proper position, as defined by synthetic spectral models, can be 
measured. Offsets in the (/ — Ha) colour appear considerably more 
frequently than those in (V — if). Therefore, we characterise off- 
sets by S(r' — Ha), the distance in this colour between the locus of 
the early-A stars in the observed photometry and their proper posi- 
tion. In the unusual case of offsets only in (/ — (') these will cause 
8(r' — Ha) to be non-zero. Observations with a measured offset of 
| S(r' — Ha) | > 0.04 were discarded, those with smaller offsets were 
corrected as necessary. The IPHAS observing includes a second set 
of exposures, displaced by 5' in both RA and Dec, taken immedi- 
ately following the original set, it was possible to verify that objects 
observed in both sets of exposures had consistent photometry. For 
all fields where both sets of exposures were returned, the photom- 
etry was found to be consistent to well within the 2% level in each 
band. 

The final sample contained observations of 132 fields and their 
offsets, in which there were 44,038 candidate early-A stars with 
r 1 ^ 19. Fig.|2]shows the coverage of the final sample, which covers 
roughly 33 sq.deg. or ~ 40% of the total area from which we have 
drawn observations. 
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Figure 3. Histograms of the ECAM distribution of the early A-stars. The 
actual observations are shown in black. The other two distributions are those 
of simulated photometry. The dashed red line indicates a model with an 
inner scale length of 2960 pc and an outer scale length of 1200 pc. The dot- 
dash green line one with an inner scale length of 2300 pc and an outer scale 
length of pc. 

able to deal with variations in metallici ty is crucial as met allicity 
affects the absolute magnitude of stars. lHales et al. (2009) found 
that for a less conservative sample of A-stars, with 8(r' — Ha) up 
to 0.07, about 15% of objects with r 1 ^ 18 were not early A-stars. 
By properly accounting for this contamination the accuracy of the 
results obtained is improved. Our simulations indicate that the con- 
tamination is, at worst, ~ 10%. The simulation includes photomet- 
ric errors and so simulated samples will contain stars with intrinsic 
colours outside the colour-cut used to select the sample, but which 
have been scatter e d in by p hotometric er r ors. 

iRobin et alj dl992h . iRuphv et all d 19961) and IRevle et all 
have all undertaken a similar approach to studying the outer 
Galactic disc. By comparing photometry from the CFHT, DENIS 
and 2MASS respectively to a Galactic model they were able to con- 
strain parameters including the thin disc scale length. However, in 
the first two cases t he observations lac k the quality and depth of 
more recent surveys. IRevle et ah I j2009t) cover a limited parameter 
space, apparently as a result of the large computational require- 
ments involved in simulating a Galactic mo del and determining ex - 
tinction using this model and the method o f lMarshall et ail 12006). 

As the reddening of a star can be determined from its position 
on the IPHAS colour-colour diagram, it is possible to correct the 
apparent magnitude of an object for this. This quantity can also be 
derived from the simulations. We refer to it as the extinction cor- 
rected apparent magnitude (ECAM) and is an estimate of the sum 
of the star's absolute magnitude and distance modulus. In this study 
ECAMs are always taken in the r' band. The difference between the 
observed and simulated r ECAM distributions is quantified using a 
% 2 test. Our aim is to find the model parameter set which minimises 
this statistic. Figure[3]demonstrates an example of poorly and better 
fitting models. 



3 METHOD 



3.1 Simulations of data 



The general approach undertaken was to compare IPHAS obser- 
vations to simulated photometry. The principle advantages of this 
method are that it allows us to take into account metallicity gradi- 
ents and low level contamination of the sample of A-stars. Being 



The Galaxy model employed i n this study uses m any of the pre- 
cepts of the Besa ncon mode l of Robin et al. ( 2003]). It is similar to 
that employed in ISale et alj J20Q9T) . In this section we provide an 
outline of the model. 
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Figure 2. The position of the 44,038 A-stars used in this study in galactic coordinates. 
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Figure 4. A graphical depiction of the type of density profile employed in 
this study. The profile is defined by: the inner scale length, L\ ; the truncation 
radius, Rj, which in this case occurs at Rj = 14 kpc, as marked by the 
dashed vertical line; the outer scale length, Lj . 



We define the disc's radial density profile with a double expo- 
nential; from the Solar Circle to some Galactocentric radius, here 
referred to as truncation radius (Rj), the density falls off with a 
scale length which we term the inner scale length (Lj). Beyond the 
truncation radius the density falls off with a different scale length, 
the outer scale length (L2). Fig[4]demonstrates the form of the den- 
sity profile we employ. We seek to measure these two scale lengths 
and the truncation radius and set constraints on the gradient of 
metallicity with respect to Galactocentric radius. 

The stellar density laws for the thin disc, in the model we em- 
ploy, are described as follows: 



P(*,*) = 



PoM)Xexp(-(0.5 2 +^ 2 /^i+z 2 A 2 ) 1/2 ) HR^Rt 
p Q /d xexp{-(0.5 2 +R 2 /L 2 + z 2 /h 2 ) 1 / 2 ) if R>R T 



Here R is the Galactocentric radius, z is the height above the Galac- 
tic plane, po the local mass density, do a nor malisation factor an d e 
is the axis ratio. Values of po are taken from lRobin et alj ( 120031) . 

The scale height, h was set at 50 pc, similar to that implied 
by the Besancon model. It is important to note that the scale height 
depends on the age of the populations studied; younger populations 
have lower scale heights. The adopted value of h is in keeping with 



measurements of young objects as follows: iBonatto et al. I d2006h 
find 48 ± 3 pc for open clusters with ages < 200 Myrs, w hich is 
an equivalent age range to that of our study; Ijoshil ( 120071) find a 
scale height of (5 7 ±4) pc for open clusters with ages < 300 Myrs; 
lElias et al.l d2006h found 34 ± 3 pc for OB stars. 

This study intentionally concentrates on the radial structure of 
the thin disc and is designed to be insensitive to the vertical struc- 
ture. Using our simulations we were able to verify that reasonable 
variations of the scale height had no significant impact on the sam- 
ple. 

In this model the central hole in the disc, which is included 
in the original Besancon model, is neglected as it has no relevance 
beyond the Solar Circle. The thick disc was included, but as it is 
old (~ 1 1 Gyrs) it has negligible influence on the results. Similarly, 
it is expected that including the stellar halo should have even less 
of a role and so it has not been included in simulations. 

We utilise the same description of warp and flare as in the 
Besancon model. Warp is described by shifting the midplane by a 
distance z W arp at a Galactocentric radius of R and a Galactic longi- 
tude of /, where: 



£warp — Zc sin(/ ^node) 
Zc ~ Ywarp ^warp ) 



(2) 
(3) 



Here / warp = 180°, specifying a node in the Anticentre direction, 
^warp = 8.4 kpc, with no warp within this radius, and Ywarp = 0. 18. 
Similarly flare is described by increasing the scale height by a fac- 
tor £fl are , which is given by: 



^flare = 1 + Yflare (R ~ ^fia: 



(4) 



Yfl are takes the value 5.4 x 10 kpc and Afl are is set at 9.5 kpc. 
As with warp, there is no flaring of the disc within this radius. 

We assume the lScalol ll986) IMF, which is valid for field stars 
in th e Galactic disc. Th e Star Formation Rate (SFR) used is taken 
(fcbm lRobin et alj J2003I) . The Galactic model assumes that neither 
the IMF nor the SFR vary with position in the thin disc. Therefore, 
these factors simply scale the number of A-stars in the simulations 
and we are not setting out to test them. This gives us the freedom 
to normalise the sample sizes of our simulated data sets as required 
and to confirm retrospectively that the implied IMF and SFR are 
reasonable. 

The mean and variation of the meta llicity of stars at the Solar 
Circle is taken from iRobin et alj J2003h . In this description, stars 
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in the solar neighbourhood with ages less than 1 Gyr have slightly 
super-solar metallicities. However, there is clear uncertainty on the 
mean metallicity of young stars in the solar neighbourhood, that 
must be borne in mind. A fuller discussion of this and its effects on 
the results we obtain is given in Section [4~4| 

The most obvious alteration is that the model is adapted to 
return magnitudes in the IPHAS r'/'Hoc system, which are un- 
supported by the Besa nco n model. The evolutionary track s of 
IPietrinferni et al.l d2004l) and lSiess, Dufour & Forestinil d2000h (for 
pre-main sequence objects) are used to convert the mass, age and 
metallicity of an object into T e ff,logg and absolute magnitude. 
These are then converted in to IPHAS colours and magnitudes us- 
ing the iMunari et al.l (2005) grid of synthetic spectra, the IPHAS 
filter profiles and atmospheric and CCD response functions for the 
INTAVFC. 

Despite the fact that ECAMs are used in this study, it is still 
necessary to include extinction in the model as the photometric er- 
ror of an object is, in part, dependent on the extinction it has suf- 
fered as well as on the distance to it. Rather than using a model 
of Gal actic extinction, such as the models of iDrimmel & Spergell 
j200lh or that used in the Besangon model, e xtinction is deter - 
mined empirically using the algor ithm MEAD dSale et alj|2009h . 
Given the extinction distance relationship found by MEAD, the 
monochromatic extinction of each star in the model is found based 
on its simulated heliocentric distances. The monochromatic extinc- 
tion is converted into extincti ons in each of the IPHAS bands using 
the R = 3. 1 extinction law of Fitzpatrick ( 2004). It should be noted 
that this conversion is not quite linear with respect to the monochro- 
matic extinction (particularly in the Ha band) and is dependent on 
the SED of the source. 

The growth of photometric errors as a function of apparent 
magnitude is taken from the observed data. The estimated photo- 
metric errors determined during pipeline processing of IPHAS data 
are fitted with a function of the following form: 

8m=A + exp (B " ,+c ) (5) 

The effect of unresolved multiplicity in star systems was in- 
cluded in the model by adding a secondary star to 56% of sys tems, 
following the binary fraction of iDuquennov & Mavoj dl99ll) . The 
addition of secondary stars was done in su ch a way tha t the final 
single object IMF was consistent with the IScalol dl986l) Galactic 
field IMF. 



3.2 Parameterisation of the simplified galactic model 

In this study we concentrated on four parameters; the gradient of 
metallicity with respect to Galactocentric radius and three govern- 
ing the stellar density profile: inner scale length (L\)\ outer scale 
length^); and the truncation radius (^7-), which is radius at which 
the scale length changes. 

Photometry was then simulated for 46,816 combinations of 
parameters, covering the following parameter ranges: 

• Inner scale length {L\): 19 values from 2300 to 3380 pc. 

• Outer scale length (L2): 16 values from to 3000pc 

• truncation radius (Rj)'- 14 values from 11 to 16.2 kpc. 

• Metallicity gradient(d[Fe/H]/rf/? c ): 11 values from -0.05 to 
-0.09 kpc" 1 . 

The model which best describes the observed data was found 
through the use of % 2 minimisation, comparing the ECAM distribu- 
tion of the observed and model photometry. For each set of model 
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Figure 5. A contour plot showing the derived confidence limits for the inner 
and outer scale lengths. Contours are at 68.3%, 95.4% and 99%. Note that 
this plot is scaled to show the entire parameter range studied and that the 
range of Z/> shown is three times larger than that of L\ . 

photometry the sample size of the observed photometry was nor- 
malised to that of the model photometry. By taking this step it is 
only the relative distribution of ECAMs which is important and not 
the total number of objects. 

The number of early A-stars at the Solar Circle is a function of 
both the local stellar density and the IMF: steeper IMFs will result 
in fewer early A-stars. By rescaling the sample sizes our results are 
not dependant on the assumed local mass densities (po) or the IMF. 

Subsequently confidence limits were calculated on the four 
parameters, by simulating 20,000 more sets of photometry for the 
best fitting parameter set and then subsequently determining which 
of the original 46,816 models fitted it best, by the same process 
of % 2 minimisation. Note that the statistical uncertainties quoted in 
the following sections for L\ ,L2,Rt and d[Fe /H] / dR are the 68.3% 
confidence intervals derived by this method. 



4 RESULTS 

4.1 Summary of the best fitting parameters 

An estimate of the inner scale length (Li) of 3020 pc with a 68.3% 
confidence range of ±120 pc. The estimate of the inner scale length 
shows little correlation with the estimates of the other parameters 
(e.g. Fig.[5). 

The metallicity gradient found is (— 0.07 ± 0.01)kpc~' . This 
is less tightly constraine d, but is well within the range of estimates 
in the literature (see e.g. lLemasle et alj2008h . 

We have been able to determine the outer scale length, L2, to 
be (1200±300) pc. The probability of the length scale being zero is 
very small, ~ 0.05%, see Fig.[6] The truncation radius is estimated 
at Rj = (13.0±0.5) kpc. The estimates of these two qualities are 
strongly covariant, see Fig. [7] so that a longer estimate of Li would 
imply a shorter ^7-. The radial density profile derived is compared 
to others from the literature in figure [8] 

4.2 Inner scale length 

As discussed in Section 13.11 our approach to dealing with stel- 
lar multiplicity is to account for it from the out set. This contrasts 
with the approaches adopted in previous studies dSiegel et aljiool: 
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Figure 6. Plot showing the distribution of the 20,000 visualisations used to 
determine the confidence intervals with respect to the outer scale length. 
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Figure 7. A contour plot showing the derived confidence limits on outer 
scale length and truncation radius. Contours are at 68.3%, 95.4% and 99%. 
Note that this plot is scaled to show the entire parameter range studied. 



Ijuric et al . 2008), where the scale length measurement assumes that 
all sources are single stars. These earlier authors preferred to ap- 
ply __an_a_p_ostoiOTi_correction to account for binarity. In particu- 
lar [siegeLetaD HqoI simulated photometry, with varying binary 
fraction, and then remeasured the scale lengths in the same man- 
ner as they had done for the real data. As the scale lengths in the 
simulations were known, the ratio between these and the lengths 
measured provided the required corrections, as a function of bi- 
nary fraction. Given a binary fraction of 50%, ISiegel et all h002h 
estimate that their measured scale lengths will be only ~ 80% of 
the true va l ue. Th ey then go on to correct their results accordingly. 
Ijuric et al.1 fc008h adopt the same correction. 

As we have included binaries in our models, our results do not 
require a correction of this form. However, to investigate the impact 
of assuming binarity, we have obtained the best fitting inner scale 
length using models without binarity. In order to focus on the inner 
scale length, we imposed a magnitude limit of r' < 16.5, so that 
only stars expected to lie within Rj are included. We found that 
the estimated values of the inner scale length using models without 
binarity were (79 ±2)% of the values found using models with 




Figure 8. The best fitting radial density profile derived in this study is shown 
in black. For comparison, the profiles of iRobin et all 1 19921) (red dashed 
line) and lRuphv etZI Jl996h (blue dot dashed line) are also shown , whils t 
the grey region indicates the inner scale length found by Juric et al. (2008). 
Note that the high Galactic latitude of the SDSS observations employed by 
Ijuric et al] feOOSD prevents them observing features far beyond the Solar 
Circle and so they cannot observe the truncation of the disc. All profiles on 
this plot have been normalised to the density at the Solar Circle. 



binarity. This fac tor is in close agreement with that determined b y 
ISiegel et alJd2002h and subsequently reused bv ljuric et aD|2008). 

If we apply this factor to o ur results, effectivel y carrying out 
the reverse of the correction of ISiegel et alj d2002h . we obtain a 
scale length of (2400 ± 80) pc, which is long er than the uncorre cted 
scale lengths of both lSiegel etall d2002h and ljuric et all d2008l) . 



4.3 Outer Scale length and truncation radius 

Our es timate of Rj = (13 .0±0.5) kpc i s substantially short er than 
that of iRobin et ail dl992l) (14 kpc) and lRuphv etaTTj 19961) (15 ± 
2 kpc). Their parameterisation of the radial stellar density profile 
was simpler than ours, in adopting complete truncation of the thin 
disc (Li = pc), as seen in Fig.[8] We find that our estimates of Rj 
and Li are anti-correlated (Fig[7J. So that if we were to assume L2 = 
pc, the best fitting distribution would be with Rj = 15.5 kpc, in 
agreement with the estimate o f lRuphvetal]dl996l). Th is su ggests 
that the simpler f ormulation assumed by IRobin et al.1 dl992T) and 
lRuphvetal.ldl996h led to the substantially longer estimate of Rj. 



4.4 Metallicity 

It has been proposed that the Galactic metallicity profile could be 
described as a step fun ction of Galactocentric radius, with a brea k 
at a radius of - 10 kpc dCorder & TwarogtoOll i lYong et al.ll2005t) . 
For our sample, such a metallicity distribution would result in there 
being a range of apparent magnitudes over which there are few 
early A-dwarfs. This is because those dwarfs immediately beyond 
the break would have discontinuously fainter absolute magnitudes 
than those in front, as a result of their lower metallicity. Such a dis- 
tribution was tested for a range of stellar density profiles and was on 
average ~ 100 times less likely than a constant metallicity gradient 
of -0.070 kpc" 1 . 

A further suggested metallicity profile i s one where th e metal - 
licity profile fl a ttens in th e outer Galaxy . ICarraro et al.1 d2007l) . 
iBragaglia et al] d2008l) and iLemasle et al] d2008l) observe such a 
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profile, whilst lChiappini et alJ ( Tl997l) and lCescutti et aljj2007h find 
it to be the natural result of their models of the chemodynamical 
evolution of a Milky Way like galaxy. In practice, the difference 
between such a metallicity profile and one with a fixed metallicity 
gradient is only manifested in small changes in absolute magnitude 
for stars at large Galactocentric radii, close to and beyond the reach 
of our observations. As a result, the method described here is not 
able to discriminate between these two profiles. 

The mean metallicity of young stars in the solar neighbour- 
hood is uncertain. Empirical age-metallicity relationships for the 
solar neighbourhood do not typically e x tend to populat i ons a s 
young as ~ 100 Myrs jHavwood 120061 : iHolmberg et al.l l2009h . 
IChenetalJj2003h measured the metallicity and distances of a cat- 
alogue of 119 open clusters. From these we make a subsample 
of nine clusters, which lie within 500 pc of the sun amd 50 pc 
of the Galactic midplane. This sample has a mean metallicity of 
[Fe/H] = -0.01. IViana Almeida etafl d2009h found that the stars 
in 11 nearby young associations have an average metallicity of 
[Fe/H] = O.OO. IViana Almeida et all [2 009) also perform a correc- 
tion to the similar data of Santos et alj J2008T) determining a mean 
value of [Fe/H] = 0.06 for a smaller sample of stars drawn from 
six local associations. 

We conservatively adopt an estimated uncertainty of 0. 1 dex 
on the mean metallicity of young stars in the solar neighbour- 
hood. The results listed above would formally indicate a smaller 
uncertainty, but our broader measure allows for systematic error 
in the determination of those values. Based on the isochrones of 
IPietrinferni et all J2004I) . a shift of 0. 1 dex in metallicity would im- 
ply an alteration of 0.06 to the r' band absolute magnitude of early 
A-dwarfs. Such an alteration would shrink or expand the distance 
to all the stars in the sample by 6% and so alter the measurements 
of L\,Lq, and Rj R . in the same proportion. Therefore, we in- 
clude this as an additional source of uncertainty on our estimates of 
these quantities, with values given in TableQ] 

4.5 Galactocentric radius of the Sun 

There is some disagreement on the distance of the Sun from the 
Galactic centre (Rq). In the study we have adopt e d the s olar Galac- 
tocentric radius of 8.07 kpc from iTrippe et alj (2008). However, 
there is substantial uncertainty on this estimate: 0.32 kpc of statis- 
tical uncertainty and 0.17 kpc of sy stematic error. In a ddition, other 
recent estimates show some scat ter: Ghez et alj j2008) found (8.0 ± 
0.6) kpc, whilst lGillessen et alj J2009I) found (8 .33 ± 0.35) kpc, by 
reanalysing the data of lEisenhauer et alj ( 120051) . Given this disper- 
sion, simulated photometry was created for a small subset of the 
original set of parameters, with Rq taking the values of 7.57, 7.82, 
8.32 and 8.57 kpc. As with the test for the effect of binarity, a mag- 
nitude limit of r < 16.5 was imposed. The values of L\ obtained 
varied by up to 2% around the estimate with the solar Galactocen- 
tric radius of 8.07 kpc. Therefore, we include this as an additional 
source of error. The effect of the uncertainty on the value of R@ on 
our estimates of the values of the inner scale length and truncation 
radius is listed in table [T] The metallicity gradient and the outer 
scale length are relatively poorly constrained, Correspondingly, the 
influence of the modest uncertainty in the solar Galactocentric ra- 
dius on these values is insignificant. 



4.6 Consistency of results with respect to Galactic longitude 

Although the values of the parameters obtained are the best fitting 
for the entire observed data set, it is not clear that this is true when 



trends with Galactic longitude are considered. For these reasons we 
split all our data (both simulations and observations) into 3 longi- 
tude subsets (160 < I < 170; 175 sC I < 185; 190 < I < 200) and 
again applied the brighter magnitude limit of r' < 16.5, to focus on 
how our estimate of L\ varies with sightline. Then, the inner scale 
length was determined with each sample. We found inner scale 
lengths of: (3030 ±320) pc, (2960 ±290) pc and (3 140 ±330) pc 
respectively for the three samples. These results are consistent with 
each other and the result for the entire sample. 



4.7 Other parameters and features influencing results 

In addition to the four parameters discussed above, there are several 
other parameters which can be varied in the model. These include 
the shape of the IMF, the form of any warp and flare of the Galac- 
tic disc and the shape of the radial abundance profile. The effect 
of varying these parameters was addressed on some of the 46,816 
para meter sets te s ted a bove. 

Re vi e et al. (2009) measure the magnitude of the Galactic 
warp to be somewhat smaller than that assumed here. However, by 
design, the form of the warp had negligible impact on the data, we 
selected the direct Anticentre for study precisely because evidenc e 
suggests it contains a node in the Galactic warp ( iRevle et al.l2009h . 

The effect of flare in the outer Galactic thin disc is to redis- 
tribute stars to greater heights above and below the midplane. This 
in turn will affect observations in two ways: the distance to the 
stars will be altered and the number which lie within a given latti- 
tude range will drop. As the intrinsic latitude range of our A-stars 
is small, the change in distance will be very small and so will have 
no significant affect on the observed ECAMs. Therefore, flare can 
only affect our results by reducing the number of objects captured 
within the latitude range we study. Using our simulations we have 
calculated that removing the flare increases the sample size by only 
~ 200 objects (~ 0.5%), almost all inside Rj. This is a result of the 
relatively low scale height of the young populations from which 
our A-stars are drawn (see SectionfSjJ. So, it is expected that only 
a dramatic increase in the flaring of the outer disc could signifi- 
cantly affect our results. For the same reason, the radial truncation 
of the thin disc that we observe cannot be the result of an incorrect 
treatment of the disc's flare. 

As discussed in section [3~7J by assumption the SFR and the 
shape of the single star IMF affects the number of stars in the 
early A range, but without altering their distribution with respect to 
Galactocentric radius. We have conducted simulations to confirm 
that the key parameters are unaffected by varying the IMF slope 
from a = —2 to —3, as expected. 



5 CONCLUSIONS 

The values of the parameters derived in this study are summarised 
in Table[T] Also listed are the uncertainties on these values, where 
the systematic component is the combination of those due to the 
uncertainty in Rp, and in the local mean metallicity. For the reasons 
given in Section l4~4l we have adopted a conservative estimate of 
the latter. As a result the uncertainties the total systematic error, to 
which it contributes, is similarly generous. 

The inner scale length of the thin dis c measured i n this 
Study is (3020± HQ statistical ± Unsystemat i c) PC Siegel et alj |2002h 



found (2500 - 3125) pc and I Juric et al] J2008h (2600 ±650) pc. 
Fig. [8] summarises som e of the different estimates obtained. Both 
Siegel et al. (2002) and ljuric et al. I d2008ri utilise samples of stars 
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Table 1. The best fitting values of the parameters, with their associated confidence intervals. 

Parameter Best fitting 68% Confidence Uncertainty Resulting From 0.1 dex Uncertainty Resulting Total Systematic 
Value Limit Uncertainty On Local Mean Metallicity From Estimate Of R@ Uncertainty 



u 
u 

R-j 

rf[Fe/H] 

Rc, 



3020 pc 
1200 pc 
13kpc 
-0.07 lqxr 



±120 pc 
±300 pc 
±0.5 kpc 
±0.01 kpc~ 



±170 pc 
±70 pc 
±0.3 kpc 



±50 pc 
±0.5 kpc 



±180 pc 
±70 pc 
±0.6 kpc 



which have large mean ages. The measurement of ljuric et all 
( 2008) is largely based on K and M dwarfs, which will have a mean 
age of several Gyrs. This contrasts with our sample of stars, which 
is predominantly young, early A-stars will only spend ~ 200 Myrs 
on the main sequence. If the Galaxy was undergoing inside-out 
formation, it would be expected that the inner scale length mea- 
sured for our sample would be longer than those obtained from the 
predominantly older samples employed in other studies. The re- 
sults obtained are consistent with this scenario, though statistical 
and systematic uncertainties prevent any direct conclusions being 
drawn at this stage. 

Additionally, we have utilised observations taken very near 
the Galactic plane, thus maximising our sensitivity to the density 
profile of the thin disc. Our estimate of L\ is longer than those of 
Robin et al. | dl992l) and lRuphvetail ( 1996), but it is unclear how 
these two studies have treated binarity. 

In common with lRobin etalj ( 1 19921) and lRuphv et al.1 dl996h 
we find that there is a knee in the radial density profile. However, 
where as they proposed that the disc is completely truncated at that 
point, we find evidence for the disc continuing, albeit with shorter 
scale length. The more gradual truncation that we observe would 
appear to be more physically reasonable, given that complete trun- 
cation would rule out both l ow level star formation out side the trun- 
cation radius (discussed in Elme green & Hunter|[ 2006) and the out- 
wards radial mig ration of stars beyo nd the truncation radius, as seen 
in the model of Rosk ar et al. (2008). 

Given our findings, the Galaxy appear s to be an examp le of a 
Type II galaxy on the scheme proposed by iFreemanl j 1970). Such 
galaxies appear t o be common in the Unive rse, accounting for 60% 
of the sample of IPohlen & Truiillcl d2006l) . We find that the ratio 
of the inner to outer scale lengths to be 2.5 ±0.6, cons i stent with 
the mean ratio of 2. 1 ± 0.5 found bv lPohlen & Truiillol d2006l) for 
classically truncated disc galaxies. 

iNaab & Ostrikerj d2006h produced a Galactic model, incorpo- 
rating many observed properties of the Milky Way. They tested the 
effect of star formation thresholds on their model and concluded 
that they would truncate the stellar surface density at a Galactocen- 
tric radius of ~ 12 kpc. This result appears to be shorter than the 
truncation radius we have measured. 

We have measured a metallicity gradient of (— 0.07 ± 
0.01)kpc~', which is w ithin the range of pr evious measurements 
in the literature (see e.g. lLemasle et alj|2008h . Frequently, the ages 
of the tracers used in previous studies have covered a wide range 
of ages, compared to our early A-stars which are between 10 and 
200 Myrs old. Therefore, our results are indicative of the metal- 
lic ity of the ISM in the re cent past, whereas other studies, such 
as lBragaglia & Tosa J20061) . have utilised older tracers, which thus 
probe the ISM at a different point in the history of the Galaxy. Our 
approach contrasts with previous studies, in that instead of using a 
small number of objects, each of which have well measured metal- 
licities, we have utilised the characteristics of a large population 
to determine a metallicity profile. We have also managed to rule 



out the possibility of th e radial metallicity profile following a step 
function (as sugested bv lYong et al.ll2005l and others). 

Fig. [9] compares observed and simulated IPHAS data for the 
same location on the sky, where the growth of photometric errors 
in the simulations have been copied from those in the observations 
and the extinction in the simulations has been derived from the ob- 
servations with MEAD. The differences between the real and sim- 
ulated observations arise for several reasons. The simulated pho- 
tometry lacks low mass (M < 0.5Mq) stars and so does not ex- 
hibit the large number of M-dwarfs visible in the observations, with 
large (r' — i ) and (/ — Ha) colours. Although the relative distri- 
bution of early A-dwarfs is not dependant on factors such as the 
IMF and SFR, the appearance of both the colour-colour and colour- 
magnitude planes is. In the future studying the distribution of stars 
drawn from a broader spectral type range should allow these factors 
to be estimated. Our simulated photometry has been derived from a 
model of a smooth axisymmetric Galaxy. On the larger scales con- 
sidered by this study such an approximation is sound, but on the 
smaller scales, such as that of an indicative IPHAS field (Fig. [9}, it 
is not. Future use of the IPHAS data will include a full mapping of 
the outer Galactic disc with all available IPHAS data over the range 
b < 5°. This will enable departures from smoothness to be clearly 
observed and lead to more focus on the vertical structure. 

This study has only utilised a small fraction of the available 
IPHAS data, with the sample limited by both spectral type and area 
on the sky. Extension to later spectral types may allow the varia- 
tion of the density profile with stellar ages to be studied, possibly 
allowing the inside-out formation of the Galaxy to be studied. It 
is currently intended to produce both observed stellar density and 
extinction maps for the entirety of this r egion. Additionall y, there 
are several other ongoing (UKIDSS-GPS lLucas et al.l2008l UVEX 
iGroot et all2 009) and upcoming (VPHAS+, VVV) optical and near 
infra-red surveys of the Galactic plane, the observations of which 
could also be used to further analyse the Galactic thin disc. 



ACKNOWLEDGEMENTS 

The authors would like to thank the referee, Annie Robin, for her 
comments that have helped improve this paper. The authors would 
also like to thank Kerttu Viironen for her helpful comments. 

This paper makes use of data obtained as part of the INT Pho- 
tometric Ha Survey of the northern Galactic Plane (IPHAS) carried 
out at the Isaac Newton Telescope (INT). The INT is operated on 
the island of La Palma by the Isaac Newton Group in the Span- 
ish Observatorio del Roque de los Muchachos of the Instituto de 
Astrofisica de Canarias. All IPHAS data are processed by the Cam- 
bridge Astronomical Survey Unit, at the Institute of Astronomy in 
Cambridge. SES is in receipt of a studentship funded by the Science 
& Technology Facilities Council of the United Kingdom. 

The simulations presented here have been prepared using the 
Imperial College High Performance Computing Service. 



10 Sale et al. 



0.8- 



0.6 

"5 

= 0.4. 
0.2 




* — — > . ■■I:...- ' — i — — t — l — 1 — 1 — 1 — 1 — 1:1 




__j 1 j j j j 1 1 j j 1_ 




■ « ... I t. ■ 1. . I 




(r' - f ) 



(r' - 1") 



Figure 9. Actual (left) and simulated (right) IPHAS colour-colour (top) and colour-magnitude (bottom) diagrams for IPHAS field 2566, with its field centre 
in the direction (l,b)=(175.6,0.2). The total area of the field is 0.28 sq.deg. The simulations lack M-dwarfs, which exhibit large (r 1 - i') and (r 1 - Hot) colours 
and more exotic objects such as Ha emitters, of which there are three prominent examples in this field. 
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